Light-mass Bragg cavity polaritons in planar quantum dot lattices 
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The exciton-polariton modes of a quantum dot lattice embedded in a planar optical cavity are 
theoretically investigated. Umklapp terms, in which an exciton interacts with many cavity modes 
differing by reciprocal lattice vectors, appear in the Hamiltonian due to the periodicity of the dot 
lattice. We focus on Bragg polariton modes obtained by tuning the exciton and the cavity modes 
into resonance at high symmetry points of the Brillouin Zone. Depending on the microcavity design 
these polaritons modes at finite in-plane momentum can be guided and can have long lifetimes. 
Moreover, their effective mass can be extremely small, of the order of 10~*mo (mo is the bare 
electron mass), and they constitute the lightest exciton-like quasi-particles in solids. 

PACS numbers: 71.36.+C, 78.20.Bh, 78.67.Hc 



I. INTRODUCTION 

The observation of vacuum Rabi oscillations between 
a quantum dot (QD) exciton and a cavity photon repre- 
sents a distinctive example of cavity quantum electrody- 
namics realized in a solid state systemi^i^i^ The strong 
light-matter coupling in networks of cavities and dots is 
very attractive for the realization of scalable quantum 
information devicesi^ Furthermore, fundamental investi- 
gations on quantum phase transitions in extended cav- 
ity systems have recently attracted a significant inter- 
esti^iii^ The strong light-matter coupling between exci- 
tons and photons is often described in terms of polari- 
ton modes, which identify propagating electromagnetic 
modes in a medium with strong dispersive properties. 
Lately, there has been a considerable theoretical and ex- 
perimental effort on polariton condensation in planar mi- 
crocavitieS] ^i^°'^^i^^ In the simplest picture, the critical 
temperature for polariton condensation is inversely pro- 
portional to the polariton effective mass. Given their 
small effective mass (about four orders of magnitude 
smaller than excitons), quantum well polaritons suggest 
the possibility of creating a condensate at high temper- 
ature. Therefore, it is important to investigate diverse 
geometrical realizations of polaritonic structures, where 
the photon and matter excitations are confined in differ- 
ent ways. These diverse geometries could lead both to 
novel quantum information devices and to novel systems 
in which coherent matter states can be obtained at room 
temperature. 

In this paper, we exploit the newly gained opportu- 
nity of growing artificial structures with spatial periodic- 
ity comparable to the wavelength of the confined quasi- 
particle, the exciton in our case, and their embedding 
in a microcavity. The main idea consists in tuning the 
exciton energy in such a way that at least two photon 
modes with different momentum are resonantly coupled 
to it due to the lattice symmetry. We assume that the 
exciton-photon coupling is stronger than the exciton in- 
homogeneous broadening. There has been many propos- 
als for engineering exciton-polariton modes using peri- 
odic structures: with periodic quantum well (QW) Bragg 



structuresfi^ii^ cavity-free three dimensional arrays of 
QDsfi^ point-dipole crystal^i^ photonic bandgaps with 
anti-dotSfii and confined QW polaritons in mesa struc- 
tures Here we consider a system with a mismatch be- 
tween the continuous symmetry of the two-dimensional 
photons in the cavity and the discrete symmetry of a 
quantum dot lattice. In the proposed structure, the po- 
lariton dispersion is entirely determined by the coupling 
between the exciton in the dots and two dimensional 
photons, which leads to in-plane effective masses which 
can be exceptionally small, of order of 10~®mo, for some 
highly symmetric points in the first Brillouin zone. The 
group velocity of 2D photons is large for large in-plane 
momenta. This can be seen as the underlying reason 
behind such a small polariton mass at the zone edge. 
The extremely small mass makes this structure a promis- 
ing candidate for novel approaches to high temperature 
polariton condensation, and for long-range coupling of 
spioi^ior exciton22, qubits. 

This paper is organized as follows: In Sec. [TI] we in- 
troduce the system and write the Hamiltonian for the 
QD lattice interacting with the two-dimensional photon 
modes. The energy dispersion of the polariton states for 
some highly symmetric points in the Brillouin zone of the 
square and hexagonal lattice is presented in Sec. HI. In 
the same section we also discuss the exciton-light mixing, 
the effective masses, and the polariton lifetime. Conclu- 
sions arc found in Sec. IV. In the Appendix we provide 
the explicit form of the exciton-photon coupling constant, 
we discuss an optimization procedure for the system de- 
sign, and we give the exact effective mass expressions. 



II. POLARITON HAMILTONIAN 

We consider the system shown in Fig. [TJ A two- 
dimensional lattice of QDs or impurities is placed at the 
center of a planar cavity structure, i.e. a Fabry-Perot 
resonator. There is currently an intense experimental 
effort for controlling the growth of QDs, based on self- 
assembly and growth on patterned substrates i^i^ Con- 
trolled impurity implantation or micro-patterned electro- 
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FIG. 1: Scheme of the system studied. An ideal quantum dot 
lattice is embedded in a planar microcavity of length Lc and 
refractive index ur. 



static traps on a quantum wel l^^'^*^ are further options to 
realize the proposed structure. We assume that the array 
consists of identical QDs and we use the effective mass, 
envelope function, and single sublevel approximationsi^ 
Furthermore, taking into account that the confinement 
potential of a QD can be modelled by a harmonic po- 
tential^ which enables the center-of-mass and relative 
motion separation;22i2^ we write the ground state heavy 
hole exciton wave function in the QD centered at Rj in 
the lattice plane as 

*j-(re, r,,) = X(R - Iij)Hp)<l)e{z,)(j)h{zh) , (I) 

where (t>e(h)ize{h}) is the electron (hole) z-confinement 
wave function, $(p) = y^2/TTa%e-P/''^ is the exciton rel- 
ative motion wave function with as being exciton Bohr 
radius, and x(R-) is the exciton center-of-mass wave func- 
tion 



(2) 



where /? is the effective dot radius. By considering only 
the ground state we assume that the energy of the first 
excited state is larger than the exciton-photon coupling, 
which is of the order of a few meV. 

We assume that the photon modes are confined by 
ideal planar mirrors in a cavity of length Lc and we con- 
sider only the lowest cavity mode with the photon dis- 
persion 
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(3) 



where c is the speed of light, is refractive index of the 
cavity, is fixed to = i^/Lc^ and q is the in-plane 
momentum. Finally, the in-plane exciton-photon or po- 
lariton Hamiltonian in the second quantization reads 



q.A 



(4) 



where hiVx is the exciton energy, crj^ (o'j,^) is the j*^ 
QD exciton creation (annihilation) operator with polar- 
ization /i, a\ q (iA,q) is the photon creation (annihilation) 
operator for a given polarization A. The constant ^ ^ 
is the exciton-photon coupling constant for the j-th QD, 
which has the property 
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(5) 



In the following we are going to focus on TE cavity 
modes, and we consider only one exciton polarization. 
For quantum well polaritons, TM modes give smaller po- 
lariton effects due to the weaker light-matter coupling)^ 
and we expect the same to be true in the lattice case. 
The lattice symmetry can be exploited by introducing 
new exciton operators 
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(6) 



After some algebra, the polariton Hamiltonian can be 
written as 



H 



= E 



qGSZ 



X]^^q+Q"q-hQ°-q+Q 
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Q 



{5q-hQ«q-hQ0-q + ^-C-} 



(7) 



where q is now restricted to the first Brillouin Zone (I''* 
BZ) due to the periodicity of the QD lattice and Q is 
a reciprocal lattice vector . The coupling constant 
expressed in terms of materials and structure parameters 
is given in Appendix VK\ 

Moreover, the light-matter interaction conserves the 
momentum q only up to a reciprocal lattice vector. The 
terms in Eq. ([7]) involving the exchange of reciprocal 
lattice vectors are known as exciton Umklapp-processcs 
in the literatureJ^i^i^ This is in contrast to the QW 
microcavity case where there is a one-to-one correspon- 
dence between cavity and exciton modes (the in-plane 
momentum is conserved exactly). They can be seen as 
interband scattering in which a valence electron with 
momentum q^^ scatters with a photon with momentum 
(\2 and creates a conduction electron Cc- If the final mo- 
mentum c\i + q2 lies outside the 1** BZ, then the mo- 
mentum of the conduction electron is flipped back to the 
first BZ by adding a reciprocal lattice vector. Since the 
photon momentum c\2 is not restricted to the V* BZ, the 
Umklapp momentum Q for excitons can have an arbi- 
trary large magnitude. This is a difference with respect 
to clcctron-phonon Umklapp processes^ where only the 
first neighbor reciprocal lattice points are involved. 

However, we remark that in Eq. ([7]) the quantum dot- 
photon coupling has a structure form factor due to the 
finite size of the quantum dots, i. e. 
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(8) 
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FIG. 2: (Color online) Scheme of photonic branches (solid 
lines) in a repeated zone scheme. The exciton energy band 
is indicated by the flat dashed line. Bragg polaritons are 
obtained by tuning the exciton energy at the zone boundary 
where photonic branches cross. Upper, Lower, and Central 
Polariton modes (indicated by UP, LP, and CP) appear due 
to the light-matter coupling. 
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FIG. 3: Reciprocal lattice of a square lattice (a) and of a 
hexagonal lattice with lattice constant a (b). High symmetry 
points discussed in the text are indicated. 



SO that we find a natural cut-ofF for the reciprocal lat- 
tice vectors. Furthermore, in the photonic case, excitonic 
corrections due to the Coulomb interaction between the 
conduction electron and the hole left in the valence band 
have to be included. 

For an exciton state with momentum q away from the 
zone boundary, many ofF-rcsonant Umklapp terms give 



corrections which do not entail qualitative novel proper- 
ties for the polariton quasi-particles. However, the situa- 
tion can be very different if the structure is built in such 
a way that the exciton is resonant with the cavity at a 
Qq which is at or near the BZ boundary. In that case we 
can choose, for instance, and so that the Bragg 
condition 



(9) 



is satisfied. In general, the number n of reciprocal lat- 
tice points satisfying this condition depends on the lat- 
tice symmetry and on the energy of the exciton. The 
mode configuration leading to the formation of Bragg 
cavity polaritons is illustrated in Fig. [2l where the exci- 
ton energy (red dashed horizontal line) is tuned in such a 
way that the polariton mixing occurs at the zone bound- 
ary. Due to the reduced symmetry of the dot lattice the 
fcz ~ 7r/ic cavity mode is folded giving rise to many pho- 
tonic branches, each characterized by a different recipro- 
cal lattice vector Q. In order to visualize this folding, the 
figure shows two repeated BZs. Bragg polariton modes 
will appear at qg where two or more photonic branches 
are crossing and are nearly degenerate with the exciton 
energy. These polaritons will have different properties 
than polaritons usually investigated at qg ^ 0. Fig. [3] 
shows the high symmetry points at the zone boundary 
for the square and hexagonal lattice for which we will 
discuss Bragg polaritons in the following. 

Due to the cut-off introduced by the finite quantum 
dot size, we can write the single polariton Hamiltonian 
at fixed q in the form 



.9q 

ffq+Qi 
ffq+Q2 

q+Q™a. 
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3q+Q. 








huj, 



q+Q„ 



7 

(10) 



Furthermore, the finite size of the dots limits the largest 
reciprocal lattice vector to |Qmaa;l 27r//3. If we keep 
only the n resonant terms at the zone boundary satisfying 
exactly the Bragg condition, the matrix Eq. pO|) can be 
diagonalized analytically and we find the two eigenvalues 



Al,2 = fu^x ± V".9qo 



(11) 



corresponding to the strongly mixed polariton states, as 
well as the (n — l)-fold degenerate eigenvalue 



A, 



(12) 



The normal mode splitting (Rabi energy), i. e. the differ- 
ence between the lowest (lower polariton) and the highest 
(upper polariton) energy at resonance, is then 



qo' 



(13) 
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FIG. 4: Bragg resonances at the M point: Ml labels the 
energy and momentum of Bragg polaritons obtained with the 
four nearest neighbors reciprocal lattice points (black dots). 
M2 is a higher order Bragg resonance involving eight second 
nearest neighbor reciprocal lattice points (crosses). 



and is proportional to the square root of the number of 
reciprocal lattice vectors involved in the Bragg condition. 

In the square lattice, we focus on Bragg polaritons ob- 
tained at the X and M high symmetry points defined in 
Fig. [3] (a). The number of photonic branches that can 
be brought into resonance is determined by the number 
of points in the reciprocal lattice equidistant to X or A/, 
according to Eq. There are two nearest equidistant 
points to X (n = 2) and four nearest equidistant points 
to M [n = 4) (see dots in Fig. [4|). However, if we do 
not concentrate on the lowest photonic branches and the 
nearest-neighbors in the reciprocal lattice, but e.g. on the 
second nearest neighbors, then a larger number of recip- 
rocal lattice vectors can be involved. For instance, there 
are eight second-nearest neighbor equidistant points to 
M [n = 8, see crosses in Fig. 01). In the following we 
will label the M polariton states as Mi if n = 4 and M2 
if n = 8. For the hexagonal lattice (see Fig. [Hb)) we 
restrict our discussion to the VF-point, which has n = 3. 
The [/-point is similar to the X-point of the square lat- 
tice and will not be discussed. Since it is possible to 
express explicitly the coupling constant for a given high 
symmetry point of length gg = ^'^'/o- (see Eq. (|B3p in 
Appendix [B| . Then the normal mode splitting can be 
written as 



nn 



(14) 



where gx is the coupling constant at the X-point. The 
values of n and ^ can then be used as scaling parameters 
to calculate the strength of polaritonic effects at different 
symmetry points, and are summarized in Tab. HI 



A. Effective mass at X-point 



TABLE I: Number of modes n and parameter ^ for the Bragg 
polaritons discussed in the text. 





X 


Ml 


M2 


W 




1 


^/2 




4/3 


n 


2 


4 


8 


3 



which is a consequence of the fact that two-dimensional 
photons at large momentum enter in the polariton forma- 
tion. In this section, we investigate this property analyt- 
ically for a Bragg polariton at the A-point. The photon 
energy in the vicinity of the X-point can be expanded to 
the second order in the in-plane momentum as 



Ex{qx,qy) 



he 

E^x + a(-27rfe/a + v^ql + ql) 



where is the A-point photon energy, a = E^/2K'^ 
and V ~ kz/K with K = \Jk'\ + (Tr/a)^. Taking only 
into account the exciton modes and the two photonic 
branches nearly resonant in the vicinity of the X-point, 
the full Hamiltonian in Eq. (|10p matrix can be simplified 
as 

gx gx 

gx a{-2'Kq.j,/a + u'^ql+ql) 
ygx a{2-Kqr,/a + v'^ql + qi^)j 

where the zero of the energy is set at E^ . The eigenvalues 
of this matrix can be found analytically by solving a cubic 
equation. From the expansion of the eigenvalues to the 
second order in q^ and qy we obtain 



E, 



UP(LP) 
^X 



± 



1 ± V2qx + 



V2 



4a:2 



{ql + v'^ql) 



f— V 

AK^gx \aKJ 



ql, 



E, 



CP 



E'x 



= 1 + 



1 



2K 



'2 



(15) 
(16) 



where the C/P, CP, and LP index identify the upper, 
central, and lower polariton branches, respectively. The 
upper (lower) sign refers to UP (LP), and 



gx 



gx 

E\- 



(17) 



is the renormalized coupling constant. From the defini- 
tion of the effective mass tensor it follows for the upper 
and lower polariton branches 



1 



E'x 
2h^K'^ 

E'x 

2n^K^ ■ 



^/2 / TT \2 

to/ 



(18) 



One of the most interesting features of Bragg polari- 
tons consists in their extremely small effective mass, 



In these expressions, we have neglected the weak depen- 
dence of the coupling constant gx on the momentum q. 
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TABLE II: Optimal values for the lattice constant a as 
discussed in App. [B] and corresponding polariton splitting 
energies h^l for several high symmetry points for two Al- 
GaAs/GaAs (1) and CdSe/ZnSe (2) systems. The cavity 
lenght Lc is 202 nm and 182 nm for the III-V and the II- 
VI systems, respectively. 



ai (nm) 
02 (nm) 
hQi (meV) 
hn2 (meV) 



This dependence gives small corrections to the effective 
mass and is discussed in Appendix [C] Since typically 
gx ~ 10~^, and Lr_ ^ a, the effective mass ratio reads 



X 


Ml 


Ah 


W 


121 


300 


383 


161 


118 


167 


373 


157 


2.2 


2.2 


1.4 


2.0 


8.8 


8.8 


5.6 


8.0 



a, 

T-yy \/2 V 



(19) 



Therefore, the polariton mass is very anisotropic at the 
AT-point. We can also compare the value of these effective 
masses with the QW polariton mass, vifhich is of the same 
order of magnitude of the cavity photon mass at q 



niph 



(20) 



and typically ~ 10^ mg. Then, we remark that for the 
the mass along the TX direction we have mxx/fTiph ^ 
10~^ and myy/niph ~ 1- This implies that a value of the 
order of 10~^too is expected. Moreover, we will see in the 
next section that isotropic masses of such a order can be 
obtained at special high symmetry points. 



III. RESULTS AND DISCUSSION 

In this section we present the numerical results for 
three high symmetry points of the square and hexagonal 
lattice. The numerical calculations include all resonant 
and off resonant terms. We have taken parameters from 
AlGaAs/GaAs and CdSe/ZnSe systems, which are typi- 
cal examples of III-V and II- VI quantum dot structures. 
In general, II-VI semiconductors microcavities have a 
stronger light-matter coupling, e.g. for CdTe microcav- 
ity the light-matter coupling is g ~ 20 meV^^ while for 
III-V systems the coupling constant is g ~ 2 meV»^ The 
cavity length, dot size, and lattice constant can be opti- 
mized in order to maximize the light-matter coupling for 
Bragg polaritons at a given high symmetry point. The 
optimization procedure is discussed in the Appendix |B] 
and the parameters are given in Tab. |TT1 

Fig. [5] (a-d) shows contour plots of the energy disper- 
sion of the upper polariton branch for Bragg polaritons 
at the J'T-point (a), M-point (6 for Mi-polariton and c for 
Af2-polariton) of the square lattice, and VK-point (d) of 
the hexagonal lattice. The energy dispersion near the res- 
onance reflects the symmetry of these points. We can also 
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FIG. 5: Contour plot of the upper polariton dispersion at the 
A-point (a), M-point (Mi in (b) and M2 in (c)), and VF-point 
(d). Note the mass anisotropy in the A-point case. 



see that the polariton effective mass is highly anisotropic 
for the A"-polariton, as discussed in the previous section, 
while the effective mass for M2 is almost isotropic. In 
Fig. [5] we show the energy dispersion of the Bragg po- 
lariton modes for X-polaritons, Mi-polaritons, and W- 
polaritons. The different behavior of the UP, CP, and LP 
modes at different points shows how richer Bragg polari- 
tons are compared to the QW polaritons. Note that the 
total number of branches crossing or anti-crossing at the 
different symmetry points is always given by n + 1 . 

At the X-point, we have n = 2 and consequently there 
are three polariton modes. In the g^-direction, the upper 
and lower AT-polaritons are similar to the QW polariton 
(see Fig. [SJd). However, there is one additional central 
mode degenerate with the photon mode along the zone 
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FIG. 6: The energy dispersion of the Bragg polaritons at selected high symmetry points. X-polaritons along VX (a) and along 
XM (b), Afi -polaritons along VM (c) and qx=qy axis (d), and W-polaritons along g^^-axis (e) and g^-axis (f). The dashed lines 
in (a) and (b) represent the uncoupled exciton and photon modes. 



boundary in the gy-direction. On the contrary, the 
dependence along VX shows the existence of (i) a lower 
branch with a negative effective mass in one direction, 
(ii) an upper branch with positive effective masses but 
high anisotropy, and (iii) a flat branch with a large mass 
(consequence of the infinite exciton mass considered in 
our model) with energy exactly between the lower and 
upper branch. 

Let us now look at the Mi-polaritons, which have five 
branches. In the qx = qy direction we can see nearly flat 
photonic branches originating from photon modes with 
Q = (27r/a,0) and Q = (0, 27r/a), while the polariton 
mixing occurs mainly with photon modes having Q = 
(0,0) and Q = (27r/a, 27r/a). A similar plot for M2 (not 
shown) gives nine polariton branches. Notice that there 
is not a saddle point at the A/-point. Polariton mode 
with isotropic negative mass (lower branch) can give rise 
to an effective negative index of refractioui^ In the W- 
point case there are four branches and the properties are 
similar to the Af-point of the square lattice. 

In the case of QW polaritons, the exitonic component 
is equally shared between lower and upper polaritons at 
the anticrossing. This is also true for the excitonic com- 
ponent in the X-point as shown in Fig. [7] (a and c). 
At the anti-crossing, the UP and LP branches consist of 
half exciton and consequently, the excitonic component 
is completely absent in the CP branch. Nevertheless, it 
increases rapidly away from the anticrossing region. This 
exciton-photon character swap is a characteristic of the 
Bragg polaritons since in all the cases investigated we 
have found analytically and numerically that the exci- 
tonic component is always equally shared between the 
lower and upper polaritons at the anticrossing, the rest 



TABLE III: Lower (LP) and upper (UP) polariton effective 
masses rrixx and rriyy at X- and M2-points for two systems 
in units of 10~^mo. For comparison, the cavity photon mass 
iTiph = hnnkz/c for the two systems are 2.1 10~^mo and 
1.7 IQ-'^mo. 



GaAs/AlGaAs CdSe/ZnSe 



X 


Ah 


X 


M2 


3.7 10"^ 


4.7 10~^ 


8.0 10"^ 


1.0 10" 


8.3 


4.7 10"^ 


6.2 


1.0 10" 


-3.7 10"^ 


-4.7 10"^ 


-8.0 10"=* 


-1.0 10" 


8.3 


-4.7 10"^ 


6.1 


-1.0 10" 



being purely photonic. 

The effective masses calculated numerically for the two 
material systems are given in Table IIIIl As expected a 
strong anisotropy at the A-point and almost isotropic 
mass at M2-point are found. The numerical values for 
the A-point are the same as calculated using the an- 
alytical formula Eq. (|18[) . thus validating our numeri- 
cal approach. Furthermore, the ratio in Eq. for 
both AlGaAs/GaAs and CdSe/ZnSe systems, for which 
gx ^ 10"'^ and^jf ~ 10"^, respectively (see Appendix [B] 
for materials parameters) is also obtained. The value of 
the effective mass in some cases is of the order of 10"^?7io 
and is extremely small for a matter like quasi-particle. As 
discussed in the Introduction, the small effective mass 
could suggest new strategies for the realization of polari- 
ton condensation in microcavities. However, the dynam- 
ics of Bragg polaritons, including the phonon bottleneck 
effects)^ and polariton-polariton scattering is expected 
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FIG. 7: The exciton component of the lower (a), central (flat) 
(b), and upper (c) Bragg polariton at the X-point. 

to be different than for quantum well polaritons, and 
needs further investigations. The small polariton effec- 
tive mass can play an important role also in quantum in- 
formation implementations. It has been recently shown 
that QW polaritons can mediate a spin-spin coupling be- 
tween two charged dots or impurities over a distance of 
the order of hundreds of nanometers due to the light mass 
of polaritouF^^ Since Bragg polariton states have effective 
masses that are two to three orders of magnitude smaller, 
the range of the spin coupling could be further extended. 

So far, the lifetime of Bragg polaritons has not been 
discussed. There is an essential difference with respect 
to the QW polariton case, where the light matter mix- 
ing occurs close to q ~ 0. The quantum well polari- 
ton lifetime depends on the photon lifetime in the cavity, 
which is determined by the reflectance of the mirrors in 
the normal direction and is of the order of picoseconds. 
However, at finite q, depending on the external index of 
refraction and on the propagation angle 9p = arccos(^), 
a cavity photon can either be transmitted through the 
mirrors, or it can remain confined due to total internal 
reflection. This means that structures can be engineered 
where the polariton mixing occurs at large q with guided 
modes, which is never possible in the quantum well case. 



Since guided photons have in principle a much longer 
lifetime, the polariton lifetime will also be much longer, 
likely in the order of micro or milliseconds, depending on 
the properties of the lateral surfaces of the cavity struc- 
tures. However, the details of the cavity modes at large 
in-plane momentum, in particular leaky modes^I due to 
the distributed Bragg mirrors, need to be carefully taken 
into account, since they will affect both the dynamics and 
the radiative properties of these polaritons. 

IV. CONCLUSIONS AND OUTLOOK 

We have investigated polariton modes in a planar 
quantum dot lattice embedded in a microcavity. As- 
suming small exciton and photon homogeneous broad- 
ening, we have considered systems in the strong cou- 
pling regime. We have studied Bragg polariton modes 
formed at high symmetry points of the Brillouin Zone. 
Truly polaritonic quasi-particles with a fifty-fifty ligh- 
matter character appear at these points in upper and 
lower branches. Analytical and numerical calculations of 
the polariton effective masses give extremely small val- 
ues, of the order of 10~®mo, and in some cases strong 
anisotropy. The small mass is due to the mixing of exci- 
tons with nearly linear two-dimensional photon modes at 
large momenta. There is a full transfer of the excitonic 
wavefunction component from the flat exciton dispersion 
to photonic branches in the strong coupling region. We 
have also discussed the possibility of polariton coupling 
with guided modes, which could lead to long polariton 
lifetime. 

Future extensions of this work could include: (i) an 
explicit theoretical treatment of the homogeneous broad- 
ening with a discussion of the strong to weak coupling 
transition, (ii) realistic structure modelling of the dis- 
tributed Bragg mirrors with leaky modes, which could 
be close to the Bragg polariton resonances, and (iii) 
sensitivity to inhomogcneous broadening effects. The 
polariton-phonons, polariton-polaritons, and polariton- 
spin dynamics for these Bragg quasi-particles also needs 
further investigation, since it is expected to be very dif- 
ferent than in the quantum well case. 

Acknowledgments 

This research was supported by the National Science 
Foimdation, Grant No. DMR-0608501. 

APPENDIX A: COUPLING CONSTANT 

The lattice exciton-photon coupling in the cavity for 
the TE mode reads^S 



= ieuj^J ^—u,,Vn^,s{0)x{<i)F{L,), (A1) 
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where V = SL^ is the volume of the cavity, S being the 
surface of the cavity, u^v is the dipole matrix element be- 
tween conduction and valence band, N is the number of 
dots in the cavity, x(q) is the center-of-mass form factor 



X(q) = x/2^/3e-i I' l\ 
and F{Lc) is the z-form factor 



FiLc) 



dz(j)e{z)(j)hiz) COS ( — 



We can assume that the QW that contains the quantum 
dots is much narrower than the cavity length Lc, which 
implies F{Lc) ~ 1. The coupling constant for TM modes 



9r 



(A2) 



is always smaller than for TE modes at finite in-plane 
momentum. We can relate the QD coupling constant in 
Eq. (jAip to the equivalent quantity in the QW case as 




QW 



TE 



qd9qw 



(A3) 



where Tiuj^'^ is a QW exciton energy, no = N/Sc is 
the QD density, af^ (ag^) is a Bohr radius of the QW 
(QD) exciton. We note that recent numerical calculations 
show that for confined excitons the effective Bohr radius 
defined by the ansatz of Eq. ([T|) is typically smaller that 
the quantum well valucj^i^ The quantity nQo2TT0^ can 
be understood as a filling or packing factor. For a fixed 
dot size and exciton energy the hexagonal lattice has the 
largest packing factor and the largest coupling at q = 
0. The value of the interband matrix element can be 
calculated conviniently from the Kane energy parameter 
using the relation^ 



2^2 



Ep 



QW 



2mn 



(A4) 



In order have a strong light matter coupling, a sys- 
tem with a large dipole matrix element and with small 
QD Bohr radius is desired. From this point of view 
II- VI semiconductors are better candidates than III-V. 
Our scheme could be also implemented with impurity- 
bound excitons, which usually have very strong oscillator 
strength, large spatial extension, and very small inhomo- 
geneous broadening4i The materials parameters used in 
the calculations are given in Tab. IIVI 



TABLE IV: Materials parameters used in the numerical cal- 
culations 



hujx (eV) 
hujS'^ (eV) 



az (nm) 



(nm) 
QQW (meV) 

TiR 

Ep (eV) 
/3 (nm) 
Ldbr (/im) 



GaAs/AlGaAs 
1.68" 
1.6 
11 
13 
1.2" 
3.5'^ 
23" 
35 
2.0 



CdSe/ZnSe 
2.50'' 
2.4 
3 
6 

2.7 

2.5'' 
21^' 
35 
2.0 



"Ref. m. 
''Ref. m. 

"Comparable with the value found in Ref. \44 
''Ref. 
"Ref. 
fRef. 



increase the polariton mixing. The coupling constant at 
a high symmetry point with go = C~ can be written as 



5,,(a,Lc,0 



1 



\/Lc + Ldbr ^ 



f3 - 
— e 



5c, 



(Bl) 



where gc is a constant which does not depend on /?, a, 
and Lc- Note that the dependence on the lattice constant 
a is non-monotonic due to the competition between the 
quantum dot density and the structure factor. Moreover, 
the exciton-photon Bragg resonance condition gives the 
constraint 



njl \ \ a J \Lc 



(B2) 



In order to take into account that distributed Bragg re- 
flectors are typically used as mirrors, we have added an 
effective Ldbr to the nominal cavity spacer length Lc 
in the coupling constant.— By expressing Lc as a func- 
tion of a from Eq. (jB2p . the coupling constant becomes 
a function of only one variable gq^, (a) . The optimal lat- 
tice spacing can then be found by analytical or numer- 
ical maximization. Defining a' = a/^, wc obtain from 
Eq. dHIl) 



gq„{a,Lc,0 ^ ^ ^9qo(o-',Lc,l). 



(B3) 



Thus one needs to optimize a only for one high symmetry 
point of each material combination and then rescale it by 
^ for other points. 



APPENDIX B: PARAMETER OPTIMIZATION 



APPENDIX C: EFFECTIVE MASS 
CORRECTIONS 



Let us now fix the QD size and ground state exciton 
energy. The lattice constant a and cavity length Lc are 
then adjustable parameters which can be optimized to 



Here we include the corrections to the effective mass 
due to the dependence of the coupling constant on q. 
This provides the exact second order expansion in and 
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Qy of Bragg polaritons and therefore the exact expression where 



for their effective mass tensor. Only 



TE 



gx close to 



the X-point is considered here. The coupling constant 
written as a function of small q around the X-point reads 



9x{qx,qy) ^ 9x 



-{27Tq^/a+ql+ql)l3^/4 



^l + {27rq,/a + ql + ql)/K^ 



Using ^ l—ia;+^a;^+0(a;'') and the expansion 

of the exponential function we obtain 



/ N Q , l.x9x , 2,x f Qx\^ , 

9x{qx,qy) = 9x+9x 7^ [jT ) +9 



2.y ( qy 



l,x 

9x 


= -9\ 


2,x 

9x 


= 9l ' 


9£ 


= 9\ ' 



TT 

2a 



1 



The exact expansion for the polaritonic modes is there- 
fore 



Si 



UP{LP) 

^x 



± 



Ecp _ -1 I 1 
E'x " 2 



,2 , 9x 47r 



9% Ka 



K 



The effective masses for the UP (LP) can be written as 



1 


El 


Tflxx 




1 


E'x 


rriyy 





9x 27r 



V2 ~9l 9'x 



(CI) 



(C2) 



while for the CP we have 
1 



I^XX 

1 



m 



yy 



E'x 

E'x 
2h'^K^ 



l.X A 

9x 47r 
.4 Ka 



(C3) 



The upper and lower polaritons small mass in the a;- 
direction is mainly due to the 0(1/50) term and the q 
dependence of the coupling constant provides only small 
corrections. The numerical calculations of the effective 
mass given in Table IIIII include the explicit dependence 
of the coupling constant on the wavevector. 
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